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Imaging the asymmetric dust shell around CI Cam with 
long baseline optical interferometry 
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ABSTRACT 

We present the first high angular resolution observation of the B[e] star/X-ray transient 
object CI Cam, performed with the two-telescope Infrared Optical Telescope Array 
(IOTA), its upgraded three-telescope version (I0TA3T) and the Palomar Testbed 
Interferometer (PTI). Visibilities and closure phases were obtained using the IONIC-3 
integrated optics beam combiner. CI Cam was observed in the near-infrared H and K 
spectral bands, wavelengths well suited to measure the size and study the geometry 
of the hot dust surrounding CI Cam. The analysis of the visibility data over an 8 year 
period from soon after the 1998 outburst to 2006 shows that the dust visibility has 
not changed over the years. The visibility data shows that CI Cam is elongated which 
confirms the disc-shape of the circumstcUar environment and totally rules out the 
hypothesis of a spherical dust shell. Closure phase measurements show direct evidence 
of asymmetries in the circumstcUar environment of CI Cam and we conclude that the 
dust surrounding CI Cam lies in an inhomogeneous disc seen at an angle. The near- 
infrared dust emission appears as an elliptical skewed Gaussian ring with a major axis 
a = 7.58 ± 0.24 mas, an axis ratio r = 0.39 ± 0.03 and a position angle 6' = 35 ± 2°. 
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1 INTRODUCTION 

The star CI Cam (MWC 84, IRAS 04156+5552), opti- 
cal counterpart of XTE J0421-I-5 60, has been inten sively 
observed since it s 199 8 outburst. iHvnes et all (|2002l ) and 
iRobinson et al.l (|2002l ) classified it as a high luminosity 
(> 10* I/q) B[e] supergiant of B0-B2 spectral type with 
an invisible compact companion (neutron star or black 
hole) and estima t e the distance to CI Cam to be >5 kpc. 
iBarsukova et al.l (120061 ) classified it as a B4III-V with a 
white dwarf companion itself surrounded by an accretion 
disc at a distance d=l. 1-1.9 kpc. The distance to CI Cam 
remains uncertain and values from the literature range from 
0.2 to 17 kpc with 5 kpc being the value generally adopted. 
The primary component of CI Cam displays all the observa- 
tional characteristics of a B[e] star: strong Balmer emission 



E-mail:nt 15@st-andrews.ac.uk 



lines, optical forbidden lines and a stron g infr a red e xcess. 
The infrared excess was first reported by I AUerJ (|l973l ) who 
found that CI Cam H-K and K-L large colour indices showed 
evidence of re-radiation from circumstcUar dust. 

The 1998 x-ray outburst of CI Cam was very unusual in 
that it was extremely fast and bright - rising to a peak lumi- 
nosity of 3 X 10"^* erg/sec (assuming 5kpc distance) within 
12 hours and decaying with an e-folding time of 0.6 days 
ISmith et al.ll 19981 : iRe^nivtsev et al.lll999l ). While this peak 
luminosity is typi cal of that reached by the clas sical black 
hole x-ray novae (|McClintock fc Remillardll2006l ). the out- 
burst is much faster than the week long rise and months 
long decay typical of these systems. On the other hand, 
the 1998 outburst was not quite as fast as those typical of 
the supergiant fas t x-ray transients (jNegueruela et aLll2006l : 
ISmith et al.ll2006l ) which contain OB supergiant stars, but 
these typically only reach 10^® erg/sec. One object which 
did have a similarly short and bright outburst is V4641 Sgr, 
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which does contain a black hole, but th e other member of 
this binary is a B9III (jOrosz et al.ll20o"ll '). not a supergiant 
B[e] as in CI Cam. The high x-ray luminosity and simulta- 
neous radio eruption have led to the suggestion that the CI 
Cam system conta inse a black hole (jRobinson et al.ll2002l : 
iBelloni et a"Llll999l ). 



CI Cam circumstellar gaseous environment has been 
subsequently probed using spectroscopy and spectropo- 
larimetry techniques. iHvnes et al.l (|2002l ) concluded from 
a series of spectroscopic observations that the dust sur- 
rounding CI Cam cannot exist in a spherical circumstellar 
shell and can only form in the outer region of its equato- 
rial gas disc. As such, the dusty environment surrounding 
CI Cam is expected to be ring shaped. From high resolu- 
tion spectroscopy observations ta ken 3 years after the out- 
burst |Miroshni^enko^^er^l] (|2002l ) showed that the emission 
lines are clearly resolved and display a triple-peaked struc- 
ture suggesting an intermediate inclination angle of the cir- 
cumstellar gas d isc which has been recently confirmed by 
lYan et al.l (|2007l ). One could expect a polarisation signature 
of the d isc orientation bu t spectropolarimetric observations 
resuhs (|lkeda et aLlbOOOD remain inconclusive. 



CI Cam's dust region, the source of the i nfrare d 
emission, has been discus sed by Miroshnichenkd (|l995h . 
iBelloni et all (|l999l ') and IClark et al.1 (120001) who mod- 
elled i ts spectral energy distribution (SEP) . iMiroshnichenkd 
(| 19951 ) derived two best models, a single BOV star with a 
1150 K dust shell and a binar y system with a BOV -I- KOII 
pair with a 1050 K dust sheU. iBeUoni et al] (|l999l ') fitted a 
Kurucz model plus thin dust model to their pre-outburst 
data and found the SED to be consistent with a ^ 22000 K 
B2 or earlier st ar plus a lar g e ma ss of hot dust. The best- 
fitting model bv lClark et"al] (I2OO0I') included a 30 00 K Ku- 
rucz model plus 1500K hot dust. Robinson et al.l (|2002l ) es- 
timated the radius of the inner edge of the dust region to 
be 13 AU< rsheii <52 AU, assuming a dust temperature 
Tdust = 1350 K and a spherical distribution of the material. 
With a distance to CI Cam of 5 kpc, this corresponds to 
an angular radius in the 2.6 to 10.4 mas range. The first 
interferometric observations of CI Cam's dust region were 
obtained with the Infrared Optical I nterferometer Array 
(IOT A) shortly after the 1998 outburst (iMillan-Gabet et al.1 
|2003| ). Their best-fitting uniform disc models gave diameters 
of 5.6 ± 0.4 mas and 7.4 ± 0.4 mas in the H and K' spectral 
bands respectively. 



2 OBSERVATIONS 

In this paper we combine near-infrared interf erometric data 
acqu i red with the free-space be am combiner (|Millan-Gabetl 
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In this paper we present the first resolved parametric 
imaging model of the hot dust region of CI Cam obtained 
using the high angular resolution provided by the IOTA and 
the PTI interferometers. We introduce our interferometry 
and photometry observations in section (2] We investigate 
the properties of the circumstellar environment of CI Cam 
through the modelling of its spectral energy distribution 
(SED) based on photometry/spectroscopy data contempo- 
rary with our interferometry data in section 13.11 We discuss 
our interferometry data modelling in section [321 And finally 
draw an updated picture of the circumstellar envelope of CI 
Cam and discuss future work in section |4l 



200 ll ) and IONIC-3 beam com- 



biner (IBerger et al. 2003 ) on the IOTA interferometer 
(|Traubl ll998l : iTraub et"Zl |2003|) wit h data obtained wit h 



the Palomar Testbed Interferometer jColavita et al.l [19991 ). 
The observations cover the period from September 1998 to 
December 2006. Figure [T] shows the (u, v)-coverage obtained 
with IOTA and PTI for the science target. The stars we used 
as calibrators are given in Table [T] 

2.1 IOTA2T free-space combiner 

These observations were carried out at IOTA using the 
two telescopes available at the time, located on the North- 
East arm. The measurements were taken in the near- 
infrared H{\o = 1.65 /xm, AA = 0.30 /xm) and K' {\o = 
2.16 fim, AX — 0.32 /jm) spectr al bands. A descrip t ion o f 
the observations can be found in iMiUan-Gabet et al.l (|2003l ) 
and iTraub et al.l (|l998h . The telescope configuration pro- 
vided a projected baseline length in the 32-36 m range. The 
starlight beam from each telescope were combined at a 50- 
50 beam splitter. The output of the beam splitter were two 
combined beams which were focused on two pixels of a NIC- 
MOSS detector. The interference fringes were produced by 
modulating the optical path difference (OPD) between the 
two beams. The beam combiner and NICMOS det ector are 
described in iMiUan-Gabet et all l|2005l,l200ll . ll999l ). 

2.2 IOTA3T IONIC-3 

Our second set of interferometric data was obtained 
using the upgrade d three-telescope IOTA interferometer 
(|Traub et al.l |2003| ) with the new int egrated optics beam 
combiner IONIC3 (|Berger et al.ll2003l ). We used three dif- 
ferent baseline configurations that were obtained by moving 
the A and C telescopes along the North-East arm and the B 
telescope along the South-East arm. CI Cam was observed 
using a standard H band filter. 

The three beams are fed to the I0NIC3 beam combiner 
which provides a pairwise combination by means of inte- 
grated optics (10). The I0NIC3 combiner contains single- 
mode waveguides. The three inputs are split with three "Y" 
junctions to provide a pairwise beam combination with an- 
other set of three coupl ers. A detailed description of the 10 
component is found in iBerger et all (|2003l ). The six light 
beams emerging from the combiner are focus ed onto 6 sepa- 
rated pixels of a PICNIC Rockwell detector jPedretti et al.l 
I2OO4I ). The optical paths are eq ualised by an automated 
fringe packet tracker described in lPedretti et al.] |200^. 

Target observations are interleaved with identical obser- 
vations of calibrator stars that are unresolved for the inter- 
ferometer. The observing p r ocedu re and data reduction are 
described in lMonnier et al] |20^[). A log of the observation 
is given in table [H 

2.3 PTI 

Our third set of data was obtained using the Palomar 
Testbed Interferometer, a three-element long-baseline in- 
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Table 1. Observation log for CI Cam. The IOTA telescope configuration refers to the location of 
the A, B, C telescopes along the NE, SE and NE arms. 



Date 






Calibrator Names 




Instrument 


Telescope configiircition 


(Adopted TJD diameter) 


1998 Sep 27 


IOTA 2T 


A35-B15 


HR1313 (1.6 ± 0.8 mas) °- 


1998 Sep 28 


IOTA 2T 


A35-B15 


HR1255 (1.2 ± 0.7 mas)*" 


1998 Nov 04 


IOTA 2T 


A35-B15 


HR1255 (1.2 ± 0.7 mas) " 


2004 Nov 26-27 


IOTA 3T 


A35-B15-C00 


HD27322 (0.24 ± 0.10 mas) " 


zU04 Uec Uz 






rliJzo94o (U.oU it U.iU mas j 








HD26553 (0.35 ± 0.10 mas) " 








HD27224 (0.87 ± 0.01 mas) 


2004 Doc 12-14 


IOTA 3T 


A28-B10-C00 


HD27224; HD26553 


2004 Dec 16 


IOTA 3T 


A28-B05-C10 


HD27224; HD26553 


2005 Nov 25 


IOTA 3T 


A35-B15-C00 


HD27224; HD26553 


2006 Oct 31 


PTI 


N-W 


HD26764 (0.5 ± 0.7 mas) " 








HD25948; HD 26553 


2006 Nov 11 


PTI 


N-W 


HD26553; HD25948; HD 26764 


2006 Nov 13 


PTI 


N-W 


HD26553; HD25948 


2006 Nov 25 


PTI 


N-W 


HD26553; HD25948; HD 26764 


2006 Dec 12 


PTI 


N-W 


HD26553; HD25948; HD 26764 



" UD diameters calculated with the fBol module of getCal (Michelson Science Center, California 
Institute of Technol ogy, |http : / / ms c.caltech.edu| 
'' UD diameter from lMerand et al.l (|2005n 
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Figure 1. The (u,v)-plane coverage of CI Cam. Squares: IOTA 
2T 1998 ; Circles: IOTA 3T 2004, Triangles: IOTA 3T 2005; Di- 
amonds: PTI 2006. 



terferometer that combines two of its three apertures to 
make measurements of fringe visibility for astronomical 
sources. For this particular observation, the measurements 
were taken with five narrow bands across the K band that 
were centred at 2.009, 2.106, 2.203, 2.299, and 2.396 ^m. For 
all of these observations, PTI's 86 m north-west baseline was 
used; details on PTI can be found in Colavita et al. (1999). 

CI Cam was observed, along with the unresolved cal- 
ibration sources HD 26553, HD 25948, and HD 26764, on 
5 nights. A log of the observation is given in Table [T] The 
PTI composite spectrometer data were calibrated using the 
WebCaliqj web interface provided by the Michelson Science 
Center. 



2.4 Photometry data 

The photometric observations used to construct CI Cam 
SED model were taken from the literature and from observa- 
tions carried out at the Michigan-Dartmouth-MIT (MDM) 
Observatory in 2004 and 2005. The data span the wave- 
length r ange from t h e UV to the infrared: iBarsukova et al.l 
(2003), IClark et all (|2000|) . Two Micron AU Sky Survey 
(2MASsT lskrutskie erail 120061). Infr ared Astronomical 
Satellite (IRAS) (|Helou fc Walkeilll988l ) and Spitzer (IRS 
spectra program ID 20702). When computing the average 
flux for each wavelength, we excluded measurements taken 
during the outburst and used post-outburst data when avail- 
able. 

The MDM photometry measurements were taken with the 
2.4 m Hiltner telescope at the MDM observatory in the UB- 
VRIJHK spectral bands in Nov-Dec 2004 and Dec 2005. 
U,B,V,R,I data were taken with the MDM 8KQ camera 
equipped with Johnson U,B.V,R, Kron-Cousins I filters. J, H 
and K observations were done with the TIFKA1V0 infrared 
camera. Table [2] summarises our measurements. 



3 MODELLING THE DUST SHELL OF CI 
CAM 

3.1 Comparison of the observed SED and 

visibihties with spherical dust shell models 

We used the radiative transfer code DUSTY developed by 
llvezic fc Elitzuij (|l997l ) in order to compute the SED for our 
models. For the sake of comparison, we constructed SED 
models for both the pre-outburst and post-outburst data. 
We computed the SEDs with the stellar radiation being 
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Table 2. UBVRIJHK photometry of CI Cam obtained by A. Tannirkulam at the MDM observatory. 



Date 
(UT) 


U 


B 


V 


R 


I J 


H 


K 


2004 Nov 27 


12.13±0.07 


12.41±0.06 


11.77±0.05 


10.79±0.05 








2004 Dec 1 










9.99±0.04 7.20±0.10 


5.68±0.10 


4.44±0.10 


2005 Dec 17 










7.01±0.05 


5.64±0.05 


4.35±0.05 



described by a black body of temperature Te//. The cen- 
tral stax is surrounded by a spherical dust shell assumed to 
begin at Ri^, inner boundary of the shell and extend up to 
Rout. The temperature at the inner boundary of the shell 
is Tdust- The dust grain type is set to amorphous carbon 
since no silicate features appear to be present neither in 
the IRAS Low Resolution Spectra (LRS) nor in the Spitzer 
Infrared Spectrograph (IRS) spectrum (D. S. Moon 2007, 
private communication). The dust grain size distribution is 
a power law with minimum and maximum size boundaries 
at a,min and a,max- The dust shell density profile follows a 
single power-law with p oc . The output parameters are: 
the angular diameter of the star 6 star, the angular inner 
radius of the dust shell Ri„, the percentage of stellar flux 
Fstar in the H and K spectral bands, the interstellar extinc- 
tion Vjb-v, the stellar luminosity hstar and the dust opacity 
in the H and K spectral bands th and tk- A rough estimate 
of the dust mass Mdust is calculated using the equation s 
which relates tv and yi-dust derived by I Hunt et al.l (|2005h . 
The values given in table [3] should be used with caution as 
the dust mass uncertainties are high du e to the inaccurate 
distance and complex dust geometry. See lHunt et al.l (|2005l ) 
for a discussion on the dust mass uncertainties. The input 
parameters and best-fitting model parameters for both pre- 
and post-outburst data are summarised in table O 

Our pre-outburst model (a) is c onstructed from aver - 
aged photometry data present ed in iBergner et al.l (Il995h . 
[12, 25, 60] IRAS fluxes (|Helou fc Walked 1 19881) "and 
IRAS LRS spectrunfl As mentioned bv lClark et"al] (20001 ) 
the IRAS spectrum is quite noisy, and was averaged before 
bein g combined with t he photometry data set. Model (a) 
uses lClark et al.l (|2000l ) best-fitting parameters as input to 
DUSTY. Their SED models, also computed with DUSTY, 
include a Te//=30000 K Kurucz model plus Td„^t = 1500 K 
hot carbon dust with a standard dust grain size distribu- 
tion and density power law index p=1.75 for both the pre- 
outburst and post-outburst SEDs. Both models differ only 
in the value of the overall optical depth of the dust enve- 
lope. The best-fitting SED model (a) is shown with the pre- 
outburst photometry data in figure [S] The model fits the 
outburst data well, however, the computed visibilities fail 
to reproduce the interferometric observations as shown in 
figure The much too small visibilities implies either that 
the size of the dust shell, Ri„ = 9.7 mas, as determined 
with the SED model fitting is larger than the size measured 
with interferometry or that the dust shell was larger prior 
to the outburst. Since no pre-outburst interferometry data 
are available we can not rule out the latter. 

We applied the same procedure to post-outburst 



photometry data. Our post-outburst models are con- 
structed from averag ed MDM (see table [2} and 2mass 
jSkrutskie et al.ll2006l ) photometry data as well as a Spitzer 
IRS spectrum. The Spitzer spectrum was downloaded from 
the Public Post-BCD Data Webserver, IRSAQ. The Spitzer 
spectrum was binned in order to obtain an homogeneous dis- 
tribution of data points across our entire wavelength range. 
Our first approach was to apply the input parameters from 
model (a) . Our best- fitting model (b) is shown in figure (2] 
The model fails to reproduce the photometry data accu- 
rately and the size of the dust shell found with DUSTY, Ri„ 
= 6.6 mas, is too large for the model visibilities to match 
our measurements as seen in figurelSJa. Our second approach 
consisted in running DUSTY for a range of T^/ / = [15000 
- 30000] K, a range of Tdust = [1000-1600] K, a power law 
parameter p between 1.5 and 2.5 and various grain size dis- 
tributions. Our motivation was to derive a model with a 
dust inner radius consistent with the size measured with 
the IOTA interferometer. A lower effective temperature and 
larger dust grain size would be expected to generate a dust 
shell with a smaller inner radius. Our best-fitting model (c) 
comprises of a T^/ / = 24000 K star and a Tdust = 1550 K 
dust shell. Our grain size distribution includes much larger 
grains up to a size of a.max ~ 5.75 with a much steeper 
density power low distribution where p = 2.35. The SED 
curve for model (c) is shown in figure [2] The fiux originat- 
ing from the star accounts for only 1.8% of the total fiux at 
A — 1.65 ^m (0.6% of the total flux at A = 2.2 pim), the 
main contribution to the total flux being the emission from 
the hot dust. The SED model is in very good agreement 
with the photometry data with the exception of the J band 
data point where the modelled flux is lower than its actual 
value. The model visibility curve still fails to reproduce the 
interferometry data with model (c) dust shell size Kin ~ 
3.3 mas (see figure [SJd). The main reason for the remaining 
discrepancy is that a spherical dust environment fails to ac- 
curately describe the geometry of the dust shell surrounding 
CI Cam as is shown in section \X2\ 



3.2 



Skewed ring model derived from the 
interferometry data 



The interferometric observations were carried out in the H 
and K spectral bands. These wavelengths are well suited to 
measure the size and study the geometry of the hot dust sur- 
rounding CI Cam. IOTA visibility data taken in 1998, 2004 
and 2005, as shown in figure |4^, agree well with each other 
which means that the structure of the dusty environment 
of the star, as seen by the interferometer, has not changed 
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Figure 2. Spectral energy distrib ution data and DU STY model of CI Cam. Left iiand panel: pre-outburst data; the model uses the 
para meters f t , T^^^t and P from lClark et ah! 1 I2OOOI ). Right hand panel: post-outburst data; the dashed line model uses the parameters 
from lciark et aL the plain line uses our new set of parameters. 




Spacial frequency [MA] Spacial frequency [MA] 

l f3k ) Pre-outburst, Spherical dust shell lf3b ) Post-outburst, Spherical dust shell 

Figure 3. The visibility curves for the DUSTY models are plotted alongside the lOTA/PTI uv-averaged visibility data points. The 
diamonds (plain lines) represent the H band data (models) while the open circles (dashed lines) represent the K band data (models). 
The left hand panel corresponds to the pre-outburst data, SED model (a) while the right hand panel shows the model visibilities for the 
post-outburst data for both model (b) and (c). 



over the years. This allows us to combine our data sets and 
use them as one in our model fitting. From the shape of the 
visibility curve shown in figure a uniform or Gaussian 
disc models are clearly incompatible with the observations. 
We measured non-zero closure phase (see figure [IJj) which 
indicates that the source brightness distribution in the H 
band is asymmetric. No closure phases measurements were 
taken in the K band. Although we display uv-averaged vis- 
ibility data in the plots throughout this paper, we used the 
original data sets in our computations. 

We tested several geometric models against the data. 
Our best-fitting model is a Gaussian skewed ring, where the 
ring emission is modulated as a function of the azimuth, plus 
a point source that accounts for the unresolved contribution 
of the cen tral star. Skewed disc models were previously dis- 
cussed bv lMonnier et aL I (|2006D in the context of the mod- 
elling of the near-infrared emission from the dust discs of 
young stellar objects. This class of models assumes that dust 



discs seen in the near-infrared will appear as a ring. The az- 
imuthal intensity modulation would either account for the 
inner edge being seen at an angle or for inhomogeneities in 
the dust distribution. If seen at an angle the ring would be 
elliptical, not circular. The parameters of the model are: the 
flux of the ring, the flux of the point source, the ring major 
axis, the ring axis ratio, the position angle of the ring, the 
skewness, and the skewed angle which is the position angle of 
the area of the ring with the greatest skewed brightness. All 
angles are given east of north. Our model uses the averaged 
H and K stellar flux derived from the SED fit to constrain 
the fiux of the two components. The best-fitting model is 
presented in figure [5] and the resulting parameters are listed 
in table g] 

We find that CI Cam hot dust can be modelled as a 
skewed elliptical Gaussian ring where the semi-major axis 
matches the inner radius of the dust found for our DUSTY 
model (c). This is quite significant since it means that the 
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Spatial frequency [MA] Position angle[°] 



||4^) Visibility data with simple Gaussian and uniform disc models Closure phase data 

Figure 4. Left panel: uv-averaged visibility data points measured with the IOTA and PTI interferometers. The solid line and dashed 
line represent the uniform disc (Syo = 5.93 it 0.09 mas) and Gaussian disc [FWHM = 3.3 ± 0.05 mas) best-fitting models respectively. 
Right panel: closure phase data measured with IOTA3T. Plotting symbols identify the separate data sets. Squares: IOTA 2T 1998 ; 
Circles: IOTA 3T 2004, Triangles: IOTA 3T 2005; Diamonds: PTI 2006 



Table 3. Input and output parameters of the SED model fitting. 
Sizes in AU and dust masses were derived assuming that CI Cam 
is at a distance d=5kpc. 



Pre-outburst Post-outburst 
Model (a) Model (b) Model (c) 



Input parameters 



Te//[K] 


30000 


30000 


24000 


Tdiist[K] 


1500 


1500 


1550 


P 


1.75 


1.75 


2.35 




0.005-0.25 


0.005-0.25 


0.005-5.75 




Output paramet 


crs 




datar [mas] 


0.026 


0.020 


0.030 


^star [AU] 


0.15 


0.10 


0.13 


Kin [mas] 


9.7 


6.6 


3.3 


Rin [AU] 


48 


33 


16 


Fstar [%], 1-65 Atni 


5.5 


1.8 


1.8 


Fstar [%], 2.2 Atm 


2.3 


0.7 


0.6 




1.16 


0.92 


0.87 


^star [Lq] 




104.9 


104.8 


TH 


0.02 


0.07 


0.38 


TK 


0.01 


0.04 


0.32 


TV 


0.07 


0.23 


0.63 


Md„^t [Mq] 


SxlO-f' 


5x10-'' 


8x10-'^ 



location of the hot dust as measured by interferometry is 
the same as the location of the hot dust as predicted with 
DUSTY for this specific model (c) which is characterised by 
a lower Tg//, a larger dust grain population and a some- 
what lower luminosity. Assuming the distance to CI Cam is 
5 kpc, we find the ring major axis to be 37.9±1.2 AU. The 
ellipticity of the hot dust region argues in favour of a disc- 
shaped dust distributions rather than a spherical dust dis- 
tribution. The skew angle is different from the position angle 
of the ring itself by about 18° which shows that the skew- 
ness doesn't come from an inner rim seen at an angle, where 
the difference in angle would be 90°, but rather arises from 



Table 4. Parameters of CI Cam Gaussian skewed ring model. 



Model parameters 



Ring flux [fraction total flux] 


0.988 


Point source flux [fraction total flux] 


0.012 


Ring major axis [mas] 


7.59±0.29 


Ring major axis [AU], d=5 kpc 


37.9±1.2 


Axis ratio 


0.39±0.03 


Ring position angle [°] 


35±2 


Skewness 


0.86±0.01 


Skewed angle [°] 


17±3 


^2 


1.1 




1.5 




1.1 



an asymmetric distribution of the hot dust around CI Cam. 
The asymmetry could arise from the formation of clumps 
during the condensation of the dust grains or, alternatively 
the asymmetric geometry of the hot dust could be generated 
by the presence of a hot undetected companion that would 
evaporate some of the dust in a non homogeneous way in 
the inner dust region of CI Cam. We also find that the ori- 
entation of the dust shell of CI Cam is similar to the ori - 
entation of the radio images l|Mioduszewski fc Rupenll2004l ) 
where they find the position angle of the radio emission to 
be about 20° east of north. 



4 SUMMARY AND FUTURE WORK 

Our SED modelling doesn't yield to any conclusion re- 
garding the spec t ral cl assification of CI Cam. According to 
lOstlie fc Carrolll (jl996h appendix E, if we assume CI Cam is 
at a distance of 5kpc, our SED modelling seems in agreement 
with a BO-BIIII classification. In order to be in agreement 
with the B0-B2 supergiant classification, the distance would 
have to be ^8kpc. On the other hand, if we assum e CI C am 
is at the distance predicted by iBarsukova et al.l (|2006l ) ie 
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U.O U.O J..U Q 

Vis2 model Closure phase model[ ] 



Figure 5. Top left panel shows the best-fitting Gaussian skewed ring model to the interferometry data. The contours in the image 
correspond to 0.1, 0.5, 0.9 of the peak intensity. Top right panels shows the model square uv-averaged visibilities (plain circles) along the 
interferometry data (open circles). Bottom left panel compares the observed uv-averaged visibility data with the model data. Bottom 
right panel compares the observed closure phases with the model closure phases. 



d=l. 1-1.9 kpc, our SED modelling is in agreement with a 
B1-B2V spectra l class ification which is in agreement with 
iBarsukova et al.l (|2006l ) Fig. 6 within the error. 

Our best-fitting model to the inte rferometry data is a 
Gau ssian skewed ring whi c h con firms iHvnes et al.l (|2002l ) 
and iMiroshnichenko et al.l (|2002l ) predictions. The skew- 
ness of the ring can originate either from the formation 
of clumps during the phase of condensation of the dust or 
from the presence of a hot companion that would evap- 
orate the dust in an non- homogeneous fashion. The ring 
is ellipt ical which confirms previous spectroscopy me asure- 
ments (IMiroshnichenko et al.l |2002| : lYan et al.l |2007| ). The 
semi-major axis matches the inner radius of the dust found 
for our DUSTY model (c) which is characterised by a lower 
Te/ / , a larger dust grain population and a somewhat lower 
luminosity with respect to previously proposed models. 

If we consider the close pair orbit in CI C am (B4III-V 
+ WD) as defined in IBarsukova et al.l (|2006l ) table 3 and 
assume that the orbit is coplanar with the dust ring, we 
can derive the semi-major axis of the orbit. From the Gaus- 
sian ring axis ratio, we find an inclination i ~ 67°. With 
asini = 48 x 10® km, we derive the orbit semi-major axis 
a = 52 X 10® km=0.35 AU=74 Rq which is weU within the 
dust ring and unresolved with IOTA3T or PTI. We don't see 



an y obvious sign of the th ird massive companion reported 
by IBarsukova et al.l |2003). The search for a possible com- 
panion in the interferometry data being complicated by the 
extended asymmetric dust emission in CI Cam. 

The ring like struct ure that we find, along with the very 
bright x-ray Fe K-line (|Orlandini et al.ll2000l ). suggest that 
the luminosity emitted by the accreting object is reflected 
by the IR ring, giving rise to the Fe-K reflection spectrum. If 
this is indeed the source of the Fe-K line, then there should 
be a delay between flux changes in the x-ray contiunuum 
and the Fe-K line, which should be equal to the light cross- 
ing time of 320 minutes if the size of the ring is indeed 40AU. 
Detection of such a delay would allow an accurate geomet- 
ric distance measure to the source by combining the x-ray 
reverberation delay and the IR interferometric angular size. 
This would be the first such distance measure to an x-ray 
binary. 

Future work should focus on the creation of a 3-D 
dust shell geometry model as well as self-consistent ra- 
diative transfer calculations. Additionally, further observa- 
tions of CI Cam usi ng the Michigan Infra-Red Combiner 
('Mo nnier et al.ll2006l) and CHARA Michigan phase-tracker 
(.Bcrger et all 120061 ') at CHARA (|ten Brummelaar et all 
l2005h will be attempted in order to acquire model indepen- 
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dent images. The high angular resolution and imaging ca- 
pability would allow us to obtain an accurate image of the 
dust region morphology with the location and dimension of 
the asymmetries detected with I0TA3T. 
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